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Abstract. Consider an ambient medium and a heterogeneous entity com- 
posed of a bidimensional material surrounded by a thin membrane. The elec- 
tromagnetic constants of these three materials are different. By analogy with 
biological cells, we call this entity a cell. We study the asymptotic behavior of 
the electric field in the transverse magnetic (TM) mode, when the thickness h 
of the membrane tends to zero. We provide a rigorous derivation of the first 
two terms of the asymptotic expansion for h tending to zero. In the membrane, 
these terms are given explicitly in local coordinates in terms of the boundary 
data and of the function /, while outside the membrane they are the solutions 
of a scalar Helmholtz equation with appropriate boundary and transmissions 
conditions given explicitly in terms of the boundary data. We prove that the 
remainder terms are of order 0(fc 3 / 2 ). In addition, if the complex dielectric 
permittivity in the membrane, denoted by z m , tends to zero faster than h, we 
give the difference between the exact solution and the above asymptotic with 
z m = 0; it is of order 0(/i 3 / 2 + |z m |). 



Introduction 

We study in this paper the behavior of the solution of Helmholtz equation in a 
bidimensional medium in transverse magnetic (TM) mode (see Balanis and Con- 
stantine [5]). The medium is made out of three materials: a central region sur- 
rounded by a thin membrane of thickness h, with 9 a curvilinear coordinate, and 
a third material, which is not assumed to be thin; see Fig. [T] This assemblage 
is submitted to a field of pulsation u>; after proper scalings, lo is included in the 
complex dielectric permittivity, which may be different in the three materials. By 
analogy with the biological cell, we call this entity a cell in an environment. In 
this article, we show that as the thickness of the membrane tends to zero, i.e as h 
tends to zero, the electric field tends to the solution of a Helmholtz equation with 
an appropriate transmission condition at the boundary between the cell and the 
ambient medium. This work is a sequel to the author's former article on the static 
case [15] . 

Let us give now precise notations. Let f2 be a bidimensional bounded domain 
composed of three subdomains: a bounded domain O c surrounded by a thin mem- 
brane Oh with small thickness h, and an exterior domain O e ^: 

n = o c uo h uo e , h . 
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Figure 1. Geometric and dielectric data. 



We suppose that the cell is strictly embbeded in the ambient domain, that is: 

(O c UO h )n0fi = 0. 

The closed curves To and Th are supposed to be of class < ^°° . We denote by To and 
respectively the boundaries of dO c and of dOh H dO Cy h- 

r = do c , 

r h = do h n do eJl . 

Let /j, e , \i m and /i c be the magnetic permittivities: they are constant and strictly 
positive. Let q e , q m and q c be three complex numbers with strictly negative imag- 
inary part and strictly positive real part; they are non dimensionalized complex 
permittivities ( see [15] or [16] for a description of the non dimensionalization). We 
define piecewise constant functions /i and q on Q as follows: 

{/He, in O e ,h, (q e ,mO eih , 
H m , in O h , q = < q m , in O h , 
fi c , in O c , [q c , in O c . 

To simplify, we denote by z the product fiq, and z e , z m and z c designate the 
restrictions of z respectively to the domains O e ^, Oh and O c . We summarize our 
hypotheses in FigQ] 
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Let 4> be a given function on dil. Its regularity will be chosen later on. We 
consider the electric field u, which solves the following Helmholtz equation with 
Neumann boundary condition: 

(la) V. ( -Vu ] + qu = 0, in D., 



(lb) d n u\ dn = 4>, on dtt. 

Here d n denotes the derivative in the direction of the vector n: n is the exterior 
normal to To, and is oriented by continuity on Th and also on dft (see Fig. [IJ. Let 
us denote by u e , u h and u c the restrictions of u respectively to the domains O et h, 
Oh and O c . These restrictions satisfy the following transmission conditions: 

(2a) — d n u c \ Fo = — d n u h \r , 

(2b) —d n u e \r h = — d n u h \ rh , 

(2c) u c |r =u & |r 01 
(2d) u e \ Th =u h \r h . 

We would like to understand the behavior for h tending to zero of the solution u 
of Problem (TTJ) . 

In our proof, we assume that fi m and q m are given constants; [i c , fi e , q c and q e 
could be continuous functions of the spatial coordinates with the imaginary part of 
q c and q e bounded away from zero, without changing the argument. 

Beretta and Francini have worked on a similar problem in [5J. They considered 
a thin dielectric material Oh in an ambient medium, and they studied a Helmholtz 
equation with Dirichlct boundary condition. They compared on the boundary of 
the domain <9f2 the exact solution to the so-called background solution defined by 
replacing the material of the membrane by the interior material. The difference 
between these two solutions restricted to the boundary d£l is then given through 
an integral involving the polarization tensor, plus some remainder terms. This 
polarization tensor is defined for instance in [3J, 0], [7|, [E], M- ^ n this paper, we 
do not use this approach since we are interested in the transmembranar potential 
(see Fear and Stuchly [11]), and in the behavior of the field in the whole domain. 
We work with bidimensional domains and we expect that the same analysis could 
be performed in higher dimensions. 

The heuristics of this work are the same as in [T5] . We parameterize Oh by local 
coordinates (r),9) varying in the /i-independent domain [0,1] x M./LZ. Here L is 
simply the length of the curve dO c 

A change of coordinates in the membrane Oh is performed, so as to parame- 
terize it by local coordinates (rj,6), which vary in a domain independently of h; 
in particular, if we denote by L the length of dO c , the variables (77, 9) belong to 
[0, 1] x M./LI1. This change of coordinates leads to /i-independent expression of the 
Laplacian in the membrane. Once the transmission conditions of the new problem 
are derived, we perform a formal asymptotic expansion of the solution of ([1]) in 
terms of h. It remains to validate this expansion. 

This paper is structured as follows. First, we suppose that the parameters q c , 
q m , q e , fi c , fi m and fi e are constant with respect to space and to h and do not vanish. 
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Moreover we assume that the imaginary parts of z c , z m and z e , which always have 
the same sign for physical reasons are negative and bounded away from 0. 

In Section [TJ we define our geometric conventions. We perform the above de- 
scribed change of variables in the membrane. We refer the reader to [15) for more 
information on the local coordinates. In Section [3l we derive formally the first two 
terms of the asymptotic expansion of the solution of our problem in terms of h. 
Section [4] contains regularity result, which is necessary for estimating the error, and 
Section [5] is devoted to estimating the error. 

In addition, in Remark[5J we give the first two terms of the asymptotic expansion 
of the electric field for a thin membrane on the boundary of C , and in Remark 1 1 01 
we consider the case z m = 0. The proofs of these asymptotics are very close to the 
proof performed in Section [5] and therefore, they are omitted. 

In the case of a biological cell, /i is identically equal to 1, z c and z e are constants 
as above and z m is very small. In Section [SJ we show that if \z m \ is small compared 
to h, we just have to replace z m by in the asymptotics found in Section [3] to 
obtain the electric field in all of the domain f2 with an error in 0(h?/ 2 + \z m \). 

Let us present now our main result. 

Main result. We give the first two terms of the asymptotic expansion of the above 
function u for h tending to zero, and we estimate rigourously the error made by 
this aproximation. 

• The th order terms. The electric fields uf. and are solution of the follow- 



(3a) 



j, aim ti 

ing problem in Q: 

{Aul + z e u c = 0, mfL\O c , 
[Au c + z c u c = Q, in O c , 

with transmission conditions 



(3b) u c \ To =Uolr > 

(3c) ^- d n u c \ To = ^- d n u e \ ro , 

and with Neumann boundary condition 
(3d) d n u^\ 9n = cj). 

• The first order terms. The fields u\ and u\ are solution of the following 
problem in f2: 

'Aul +z e u\ = 0, in fl\O c , 
(4a) I Aul + z cu\ = 0, in O c , 

with the transmission conditions 

(4b) — <9„ui|r - — d n ul\r = ( — ) <9 t 2 Uo|r + {q m - qe)u L Q \r , 

(4c) ullr, - u{\r = Me — d„<u%\r . 
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We have the following theorem. 



Theorem 1. Let O c be a bouded domain with smooth boundary, and denote by k 
the curvature of dO in local coordinates. Let ho be such that: 

1 



< h < 

\\K\\oo 

Let h belong to (0, h ) and 4> to H s (dtt), s > 7/2. 

We denote by u the solution of Problem fTJ). Define (v e ,v c ) in fl by: 

v e = Uq + hu\, in£l \ O c , 
v c = Uq + hul, in O c . 

Then there exists an h-independent constant C > such that 

IK-w c |Ui( 0e )<C/i 3/2 ||^IU S (an), 



\u-v e \\ L 2 [0h) + 



— Vu - — Vv ( 



<Ch 3 / 2 \\ 



L 2 (O h 



\\u-v e \\ m(0eh) <C/i 3 / 2 ||0|| ff3(af2) . 



1. Geometry of the problem 

The boundary of the domain O c is assumed to be smooth. The boundary To 
is counterclockwise oriented, and we denote by dt the tangential derivative along 
O c . Thanks to a change of units of length, we may suppose that the length of To 
is equal to 2ir. We denote by T the flat torus: 

T = R/2ttZ. 

Since To is of class , we can parameterize it by a smooth function 'J from T to 
R 2 satisfying: 

V6»eT, |*' (0)| = 1. 
The following identities hold: 

r o = {*(0),0eT}, 

and 



Tfi = {^(9) + hn(8),6 e T}. 

Here n(8) is the unitary exterior normal at \P(0) to IV We parameterize the 
membrane Oh as follows: 

O h = {#(0) + hnn(9), (n, 8) e]0, l[xT}. 

We define now: 

$(77,6!) = ^(9) + hnn(9). 

We denote by k the curvature of the To in curvilinear coordinate, and by 8. the 
curvature in Euclidean coordinates: 

(5) VxedO, £= Kof^i). 
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Let ho satisfy: 

(6) < h < F - ir 



Then, for all h in [0, ho], there exists an open interval / containing (0,1) such 
that $ is a diffeomorphism of class < rf°° from I x K/27rZ to its image, which is a 
neighborhood of the membrane. The metric in Oh is given by: 

(7) h 2 d v 2 + {l + hr 1 K(6)) 2 de 2 . 

We use two systems of coordinates, depending on the domains O e h, O c and Oh'- in 
the interior and exterior domains O e ^h and Cl we use Euclidean coordinates (x, y) 
and in the membrane Oh, we use local coordinates with metric (J7J. 

2. Statement of the problem 

In this section, we express our Problem (TTJ) in local coordinates. It is convenient 
to write: 

V6>eT, $ O (0) = $(O,0), $i(0)=#(l,0). 

Let us denote by u e and u c respectively the electric field in O e .h and in O c , written 
in Euclidean coordinates, and by u m the electric field in Oh in the local coordinates: 

u e = u, in O e< h, 
u c = u, in O c , 
u m = 1/0$, in [0,1] x T. 

We have shown in |TS] that the Laplacian in the local coordinates is given by: 

Therefore, we rewrite Problem |T|) as follows: 
(9a) Au e + z e u e = 0, in O e , h , 
(9b) Au c + z c u c = 0, in O c , 
V(i),S)e[0,l]xT, 

(9c) d v H-^Hld^A + d e ( y Y^- K deum ) + Zmhf{1 + ^ K )" m = °' 
with transmission conditions (J2|) expressed in local coordinates at r\ = 0: 



(9d) — 9„m c o$ = TT-ftju" 
(9e) u c o$ = u m |, =0) 
at T) = 1: 

(9f) — <9„u e o$i = t—W 
(9g) « e o$i = u w '| 1J=1 , 



»7=0 



»7=1 
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and with boundary condition 
(9h) d n u e \ dn =0. 



3. Formal asymptotic expansion 

In this section, we derive asymptotic expansions of the electric field (u e , u c , u m ) 
solution of (|9|) in terms of the parameter ft. In the limit, we want to be able to 
replace the membrane by transmission conditions. 

We multiply (I9cp by h(l + hr\n) 2 and we order the result in powers of h, in order 
to obtain the partial differential equation (PDE) satisfied by u m : 



V(r?,0)e[O,l]xT, 







2„,m 



U 



tin {irid^u 711 + d n u m ) + h 2 {3r) 2 K 2 d*u r ' 
2T 1 n 2 d ri u m + 8 2 u m + z m w m } 



(10) 

+ /i 3 {?7 3 K 3 9 2 u m + f] 2 K 3 d v u m + 7]Kd 2 u m 

- 7]K'd g u m + 3z m riKu m } + 3h 4 z m rj 2 K 2 u m + h 5 rj 3 K 3 z m u m = 
We assume the following ansatz: 
(11a) u e = Uq + hul + ■ • • , 

(lib) u c = Uq + hu\ H > 

(11c) u m = ug* + ft< + ■■■ . 

We will to derive the first two terms of the asymptotic expansions of u e , u c and 
u m by identifying the terms of coefficients of a given power of ft. 

We extend formally u e to tt\O c , by extending a finite number of coefficients of 
the powers of ft. Moreover, we suppose that <p is as regular as needed. We will also 
need the first two terms of 

u e o $(77, 6>) = u e o (tf (0) + hr)n(0)) 

and d n u e o $(77, 0). This amounts to composing two asymptotics series. We remem- 
ber that we introduced $o = ^ an d $i=<I>(l,-)to homogeneize our notations. A 
simple calculation gives: 

u e o $ = u e o $o + ft (u\ o $o + r\d n u% o <& ) H , 

and similarly 

d n u e o $ = <9„u e o $ + h (d n u\ o $ + V 9 n u o o $ ) H • 

These expansions enable us to rewrite transmission condition ([9f| as: 

^ ( a„ug o $ + h {d n u\ o $ + flgug o $ ) + • • • ) 
(12a) Me V y 

= Win + ft + ft 2 fl,Hj m | H + • • • , 

and transmission condition ( |9g| ) as 

(12b) UoO$ + ft Ko$ + d n u e o<P ) + ■■■ = + fturU=i + ' ' ' • 
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Observe that we have chosen to limit the order of explicit asymptotic expansions 
to what will be needed below. 

We systematically substitute the fields u e , u c and u m by their asymptotic ex- 
pansion (jlip in ([9]). For transmission condition, at 77 = 1, it is more convenient to 
use transmission conditions (| 1 2|) instead of (|9f[) -(|9g|). 

We are going to select all terms of an appropriate order in these expanded equa- 
tions in order to get the conditions satisfied by it™, uf, u\ and u 1 ?\_ 1 (i = 0, 1). 

First step : identification of O*' 1 order terms. Substituting into (TT0|) the field 
u m by its expansion (|llc|l we obtain: 

(13) d^u'o — 0, V(?7, 9) e (0, 1) x T. 

Moreover, we obtain easily: 

(14a) Au e Q + z e u e Q = 0, in fi \ W c , 

(14b) Au c + z c u c = 0, in O c , 

and the boundary condition 

(14c) d n Uo\ 9 n = <f>. 

Equality (|9dj) implies: 



d v u™\ v=0 = 0, 

and equality (|12ap implies: 

d v u™\ v=1 = 0. 

Therefore, depends only on 9. By identifying th order term in (|9d |) ~ (|12ap . we 
infer: 

(15) <o$ = = <o$ , 
thus we obtain the following transmission condidtion: 

(16) wgo$o = w§o$ - 
We will determine later on. 

Second step : identification of first order terms. Substituting into equal- 
ity (|10|) the field u m by its expansion (|1 lcjl . and using that in (0, 1) x T, we have 

we obtain: 

(17) d 2 v u? = 0. 
Moreover, we obtain easily: 

(18a) Au\ + z e u'l =0, in n \ ~0~ c , 

(18b) Au\ + z c u\ = 0, in O c , 

and the boundary condition 

(18c) d n ul\ an = 0. 

Equality ([9d]) implies: 

(19a) fy«iVo = — s »«o°*o, 
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and equality (|12ap implies: 

(19b) ^1^=1 = 0*0. 

Me 

We infer the following transmission condition between (9„ug and <9„ug: 
(20) — itg o $ = —«o $ o- 

Mc Me 

Therefore, with flip. ([TO]) and we infer that (ug, ug) satisfies the following PDE 
in f2: 



(21a) 



Awg + z e ug = 0, in O e j t , 
Auq + z c Uq = 0, in O c , 



with the transmission conditions 

(21b) u c \r = ug|r , 

(21c) — d n Uo\ ro = — a„Uolr » 

and with Neumann boundary condition: 

(21d) d n u e 1 \gn = 0. 

According to (fT5]l . is equal to: 

(22) VM)e[0,L]xT, t#(»,,d)=u8o#o(0). 

We have determined ug, ug and Uq". 

Observe that the identification of the first order term in (|9e|) implies: 

(23) ur| n=0 =^o$o, 
and will be determine later. 

Third step : identification of second order terms. According to (fT7|) and 
([22]) . we have 

= eg< = 0. 

Therefore, by identifying the second order term in h of (110|) . we obtain: 

(24) 9^ u m + mi = Qj 

with 

(25) m x = Kd v uT + d 2 e u™ + z m v%. 

Observe that mi depends only on 9. 

The identification of second order terms of (|9d|) implies 



and those of (|12a[) implies 



d v u™\ v =o = —d„u c 1 o<i>( ) , 



(26) $,«aVi = — o $ + 9>g o $„) 

Me 
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Observe that <9,,u™ depends only on 9 thanks to (fTT]) . Thereby integrating ([24j) 
with respect to 77 we obtain: 

(27) fyv? = -Wii + — 3„^o$ . 

Me 

From (|26p and (|27p we will obtain a transmission condition for (u\,u\). More 
precisely, taking r\ = 1 in ([2"7| , with the help of ([2T>]) we obtain: 

(28) — d n u\o<S>Q - —d n u\o<S>v = m x + n m — 9X°*o, 

Me Me Me 

and the right-hand side of (|28| is entirely determined. 
By identifying the terms of order 1 in (|T2b"|) we obtain: 

(29) u;o$o-«fo*o=fl-^)fiU8o$o- 

For convenience we write equations satisfied by (u\,u\) in Euclidean coordinates. 
Thanks to fTI|). Ip8| and ([211) we mfer tnat ( u i> u i) solves: 

{Au^ + z e uf = 0, in Q.\Ol, 
Au\ + z c u\ = 0, in O c , 
9 nUf\ an = 0, 

with transmission conditions 

Mm ^ c 1 Mm o e 1 a2 c 1 el 

— 9 n u 1 \r 9 n u 1 \r = d t u \r + z TO Wo|r 

(30b) Me Me 

+ — «|r + — ^,u§|r , 

Me Mc 

(30c) «f|ro-«!|ro = (l-^)ft,«5|ro- 

In Section |H we prove the existence and uniqueness of {u\,u\) defined by ([30]) . 
Remark that u™ given by equality (j23j) is entirely determined and c^it™ is entirely 
determined by (|27l) . 

Observe that (j30bj) contains a second normal derivative; this is a feature of the 
asymptotics of a cell in an ambient medium; no second derivative appeared in [15) , 
where there is a cell with boundary condition on the exterior of the membrane. Let 
us summarize the first two terms of the asymptotics we obtained formally. 

• The th order terms. The electric fields Uq and are solution of the follow- 
ing problem in Q: 



(31a) 



Au e + z e u c = 0, in ft\O c , 
Auq + z c u% = 0, in O c , 



with transmission conditions 



(31b) uo|r = "olr , 

(31c) — <9„u§| ro = — 8„Uq\ 

f^c Me 

and with Neumann boundary condition 
(3 Id) d n u^\ 9n = (f). 
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In the membrane, the field u™ is equal to: 

(32) V(t7, 9) e [0, 1] x T, u% = ugo* o (0). 

• The first order terms. The fields u\ and u\ are solution of the following 
problem in Q,: 

!Auf + z e u\ = 0, in O \ O c , 
Auf + z c wf = 0, in O c , 

with the transmission conditions 
1 „ „, 1 „ „, 1 



(33b) 



<9n<|r - — d n u\\r = — (<9 t 2 <lr + z m u c Q \r Q ) 
+ — S£t*olr + — -fi^ttglro, 



(33c) uc^.^^^j^^i 



In the membrane, we have: 



(34) V(J7, 0) e [0, 1] x T, u? = T)—d n ul o $ + it? o *„. 

Mc 

Remark 2. FTe maj/ write 9^ug|r in terms of d n UQ\r , o/ug|r and of its tangen- 
tial derivatives. Actually, we perform the change in local coordinates in a neighbor- 
hood of dO c . According to (|31|) . the following identity holds along Tq: 

dn u o\r a = -Ad n u e \ To - <9 t Xh - z e u e \ Toi 

thus we may rewrite transmission conditions (|33bj) (|31c[) as follows: 

— 9 n Ui\r - ^-d n ul\r = (— ) d?u c Q \r + (q m - q e )u L Q \r , 

Mc Me \ pm Me / 

hence l|3"Tj) - (j3"3"|) are equivalent to ©-([l]). 

We have given the first two terms of the asymptotic expansion of u e 7 u c and u m . 
It remains to prove that the remainder terms are small. First we need to study the 
regularity of Uq and u\ in a neighborhood of T . 



4. Regularity Result 

In this section, we study the regularity of the solution of Hclmhotz equation 
with our transmission condition, which is non usual. This result is required to 
prove Theorem Q] of Section [5l which estimates the errors between the asymptotics 
and the exact solution. The following result is natural and expected; it is very close 
to a result of [12] (Appendix, page 147) by Li and Vogelius, but different enough 
to require a proof. We thank very warmly Michael Vogelius for his suggestions on 
the reflection principle. 
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Theorem 3. Let G belong to H S (T ), s > -1/2. Let (U e ,U c ) be the solution of 
the following problem: 

V. (^- VC/C ) + 1cU c = 0, in O c , 
V. (j-S7U e ^j + q e U e = 0, in n \ O c , 

with the following transmission condition: 
U e \r B =U c \ ro , 

—d n U e \ r<) -—d n U c \r =G, 

and with the Neumann boundary condition on dQ 
d n U e \ dn = 0. 

Then we have: 

U e €H s+3 / 2 (Q,\0~ c ), U c g H s+3/2 {O c ). 
Moreover let m be a non negative integer, and s > m + 1/2. Then, 



U e £ c if m (n \ O c ), U c £ tf m (O c ). 

Proof. Since To is smooth, we use local coordinates in a neighborhood of IV Ac- 
tually, as in Section [T] there exists hi such that: 

n = {¥(0) + h lV n(9), (77, 6) e (-1, 1) x T} , 

is an open neighborhood of Tq and 

(rj,6) ^ y(9) + hir)n(9) 

is a diffeomorphism from (—1, 1) x T to We denote by g the function G written 
in local coordinates: 

V6»eT, g(6)=Go^{6). 

We denote by C the cylinder [0, 1] x T and by (C) the space of the functions 
a defined on C such that: 

1 P 2ir 

2 



\ 1/2 



1 r2 * / l + h 1 riK . m|2 hi 



Jo 



\d n a(v,0)\ z + —^\d e a(r ] ,9)\ 2 )dr ] de 

til 1 + hiTJK 



is finite. We equip H} n {C) with such a norm, which is equivalent to the ordianry 
norm 

12 1/2 

f * (H 2 + |d„a| 2 + \dga\ 2 ) d9dn) , 
'o Jo J 

because we have the following equality: 
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We use a partition of unity and classical elliptic regularity to reduce our problem 
to establishing the regularity of the solutions (V e , V c ) of the following problem: 

V(»/,0)e[-i,o]xT, 
n + hivn \ + ds f hi d \ + + h = ^ 

V h iHc J VI 1 + hir)K)(i c J 
V(?7, 0) e [0, 1] X T, 

+ \ + f hi dgye \ + x + = 0j 

V mfl e J \{1 + hiT]K)fl e J 

with Dirichlct boundary conditions 

F c |„=-i = 0, V e \ v=1 = 
and with transmission conditions 

V c |,, = o = V e \ v =o, 

—d v V e \, 1=0 - —d n V c \ n= o = g, 

We use the reflection principle, suggested by Vogelius and coming from an idea 
of Nirenberg (see [12], page 147 or |T] and [2]). With the help of this principle, 
we transform transmission conditions into boundary conditions. We define V r on 
[0, 1] x T by: 

VM)€[0,1]XT, V r (r ] ,d) = V c (-r ] ,6). 
The functions V e , V r satisfy the following problem in (0, 1) x T: 
V(r?,0) e (0,1) x T, 

(35a) d v ( 1+ u hir]K d v V e ) + do ( fel d e V e ) + q^l + h lV K)V e = 0, 
(35b) d v ( 1 ~ hlVK d v V r ) + dg (- fel d e V r ) + q c hi(l - h lW )V r = 0, 

with Dirichlet boundary conditions in r\ = 1 

(35C) yVl=0, V e \r,=! = 0, 

with boundary conditions in r\ = 0: 

(35d) V t \ v =q - V e 1^=0 = 0, 

(35e) — a,V|^ + — d n V% =0 = g, 
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Multiplying (|35ap by V e and (|35b[) by V r , integrating by parts and summing, we 
obtain: 



^ \d V e \ 2 + - !^-^\deV r \ 2 - q c hi(l - h lVK )\V r \ 2 



(36a) 



{l + hirjn)^ e (l-hirjn)fj, 

q e h x {l + h lVK )\V e \ 2 \d V d9 = I ( -^^y e |„ = oF"^=o 



2tt / ^ 



\hlfi e 



+ -^—d v v r \ v=0 v r \ v=() 

tliflc 



Using boundary conditions (|35d[) - (|35ep , we obtain: 



(36b) / — d v V e \ v=0 T\ v=a + — d n V r \ v =oV Vo = / 5^ e d0. 

JO VMe / JO 

We argue as in [14] or in [4] , and the reader will verify that (136a[) and (|36bj) suffice 
to give existence and uniqueness of solutions of (j35j) in iJ^(C). 

To obtain the regularity result, we just have to apply the method of frozen 
coefficients. Let 8q S T, and denote by kq the value of k at ^o- A classical argument 
(see for instance pQ, [2], [10], [13] or [14] ) shows that (F e , V r ) have the same 
respective regularity as (V , V") solution of: 

VM)e[0,l]xT, 

All 1 + /ll?7Ko (1 + ftll/Ko) 

L&v" ^ a v" + hld2eV " = o, 

h\ n 1 — hir]Ko 71 (1 — hir]Ko) 2 
with Dirichlet boundary conditions 

vVi = o, y"|,=i = o, 

with transmission conditions 

V |»j=o = V \ v =o, 

—d v V'\ v=Q + —d v V"\ v=a =g. 

Me f^c 

The regularity results of this last problem is obtained directly by working in Fourier 
coefficients, hence the regularity result, in ~f\. The end of the proof follows by 
classical regularity theorems (see [14] for instance). □ 



5. Error Estimates 

We give an error estimate, which proves that the first two terms obtained in 
Section [3] through a formal argument are indeed the first terms, i.e. the remainder 
is smaller. 
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Remark 4. Recall that the L 2 not 
denoted by \\cx\[\ai,% (c), is equal to: 



of a 0-form a in C with the metric 0, 



Afi^(C) 



2tt 



h(l + hn K )\a(r],e)\ 2 drjd9, 



o Jo 



and £/iai tte L 2 norm of the exterior derivative da of a, denoted by ||da||AiL 2 (C) 
is equal to 



I da 1 1 a 1 14(c) 



1 + hrjK 2 
7 \d v a{r),6)\ 



1 + hrjK 



|9 e a(7 ? ,0)| z dr? d0, 



IV (aor 1 ) || 2 L2(0h 



Let us prove now TheoremQ] We remember that u is the solution to Problem ((T|) 
and that (v e , v c ) are defined in Q by: 

(37a) v e = ul + hu\, in Q\O c , 

(37b) v c = ug + Zmf, in O c . 

We have to prove that there exists an /i-independent constant C > such that 
(38a) \\u-v c \\ m{0c) < Ch 3 / 2 U\\ H s (dn) , 



(38b) w e || A o L ^ (c) + 

(38c) 



du dv 



< chV 2 w\\ H . m , 

\\u-v e \\ H i (0i! h) < Ch 3 ' 2 \\(t>\\ H s {m) . 



Remark 5. The estimates of Theorem Q] are piecewise H 1 estimates since esti- 
mate ()38b[) involves 

—du —dv e , 

f-m Me A 1 i=,(C) 

which is not the norm of a difference of gradients. However, we could have global 
estimate with an appropriate norm involving the permeabilities [i c , \i m and \i e and 
by defining an appropriately modified fi; details are left to the reader. 

Since <j> belongs to H s (dQ), s > 7/2 and Theorem H holds, u e Q belongs to "if 3 (Oh) 
and u\ to ^ ,2 (C; l ). To prove Theorem [1] we need the following lemma. 

Lemma 6. Let h belong to (0,ho)- 

Let Uq, u\, u™ andu™ be defined by pip and (|33[) . We denote by v the following 
function: 

(39) V(ry, 6) £ [0, 1] x T, v{q, 6) = u^ + huf. 
Then, there exists an h-independent constant C > such that 

(40) \\v e o^~v\\ A o L 2 {c) <Ch^ 2 U\\ H s {anh 



(41) 



— d (v e o$) — dv 



A!i 2 (C) 



< ch^ 2 u\\ H s [dn) 
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(42) 



||u e o$i - v\n=i\\A°L2(j) < Ch 2 \\(j)\\ HHan) , 
\\d g v e o<$> 1 - dev\^i\\ A o L 2 m < Ch 2 \\cj)\\ HS{dn) , 
. ||a|w e o$i - dgV\ v= l\\ A L 2 (T) < Ch 2 U\\ H s (dn) . 



Remark 7. A priori, one would have expected that the statement of Theorem [7] 
would have given a comparison of the exact solution with its asymptotics in the three 
regions O e ju Oh and O c . Actually, Lemma\^ shows that we may dispense with the 
asymptotics in Oh, provided that v e defined by (|37|) has been extended up to the 
inner boundary of the membrane, and this is precisely how v e has been constructed. 

Proof of Lemma\Q Since UqO& belongs to ^ 3 ([0, 1] x T) and since u\ o$ belongs 
to ^ 2 ([0, 1] x T), using Taylor formula with integral remainder, we obtain for all 
M)e [0,1] xT: 

v e o$| (t7i9) = Uo 0$ l(o,e) + hr/d n u e o<£\ i0 . g) +hu\o^\ {0j6) 

+ h 2 V 2 J (1 - t) (d n u\ o 9) + ^-^d 2 n ul o $(pn, 9)\ dt, 

and 

a, ; (w e o$)| (i)e) = /id n it§o $| (0 , 9) +h 2 (^d n ulo<S>\ {ve) 
Since we have: 

vfa, 6)=u c Q o $ (9) + hrj^d n u c o $ o (0) + hu\ o $ (9), 
d v v(n,9) = h^d n u c o$ (9), 
using transmission condition (|31b[) , we obtain for all (rj, 9) £ [0, 1] x T: 

(v e o $ - 5)1,^ = hL(l-^j d n u% o $ + «f o $o - «i o $0 

+ to? jf V - t) ^„uf o 9) + r?^^^ o 0)) dtj . 

This equality implies directly estimate (|40|) . Moreover, using transmission condi- 
tions (|31bp and (|33c|) . we obtain for all 9 £ T: 

v e o $! (0) - 57(1, 9) = h 2 J (1 - t) (d n ul o *(t, 0) + fc^^ug o $(t, 0)^ di, 

which implies the first estimate of (|42|) . Applying the same reasoning to <9#i; e o<I> 
and to dgV e o&, we obtain the two last estimates of (|42|) . 
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Observe that: 

—dr,(v e o^)(r,,0)- —d,v(v,0) = —(O)lvf [ (l-t)dlu e o d< 

Me Mm Me \ JO 

hence estimate (|4Tj) . This ends the proof of Lemma [6] □ 

Let us prove now Theorem [TJ 
Proof of Theorem [7J Define mo by 

Mm n c t 
m = d„u o$ , 

Me 



and u™ by 



?? 2 



where mi is defined by (|25[) . 
Let 

(43a) VK e = u e ~ (u% + hu\), in O e ,h, 

(43b) iy c = u c - («g + ftuf) - h 2 5 c , in O c , 

(43c) VF m = u m — « + /zuf) - /i 2 B m , in [0, 1] x T, 

where 

(44) B m ( V , 6) = uf + a m (9) + r}b m {rj), 

and B c , a m and b m are allowed to depend on h and will be chosen later, so that 
yield the easiest estimates of W e , W c and W m . 

Let us write the problem satisfied by (W e , W c , W m ). In order to simplify the 
notations, we introduce Jzf, the Hclmholtz operator written in the local coordinates 

(V, ®) S iven b Y 



1 + 

ft 



We obtain 



(45a) AW e + z e W e = 0, in O e , h , 

(45b) AT^ C + z c W c = ~h 2 (AB C + z c B c ) , in O c , 

V(r),6)£ [0,1] xT, 
iff" 1 = -J5? (ttg* + /i< + ft 2 B m ) , 
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with transmission conditions coming from |2j) 

—d n W c o^>o d ri W m \_ +h(—d v B m \ ri =o-—dnU c 1 o$o 

(45c) Me ft/im ' Vm Me 

-h 2 d n B c o$ , 

(45d) ^ c o$ = W"|„ =0 + & a (#"Vo - S c o$ ) , 

(45e) —d n W e o$i =-J—fd v W m \ + hd v u? + h 2 d v B r 

(45f) - — 9 n ugo$i - — a n «fo$i, 

Me Me 



1 77=1 



(45g) 



^ e O$! = ^ m |, )=1 +/ l 2 B m | I)=1 +<Vl 
+ /lU™l»7=l - Uo ^! ~ /lUi 0$1, 



and the boundary condition 
(45h) d n W e \ m =0. 

We calculate JfW m , knowing that d^u™ and d^u™ vanish and we obtain: 

^ SfW m = - hnd, lU r i - hd n ((1 + hrjn^B" 1 ) - hdg [ j^h^, de J ( u ™ 

+ hu 7 ? + h 2 B m ) - z m h (1 + hrjK) (tig 1 + hu? + h 2 B m ) , 
and we find that in the above expression, the coefficient of terms of order 1 in h is: 
(47) -d 2 B m - Kd v u™ - d 2 e v% - z m u™. 

By definition (|4~4"|) of B m , the first term of (|47[) is mi, and according to the defini- 
tion ([25]) of mi, the expression (|47|) vanishes. 

We will determine a m and 6 m so as to have nice transmission conditions. Observe 
that if y satisfies, in the weak sense: 

^7yj +qy = 0, mO e>h UO h , 

with discontinuous fj, on the outer boundary of the membrane Oh , then the trans- 
mission conditions on this boundary are: 

y e o<5> 1 = y m \ v=1 , 

Me '~ 

Therefore, it is natural to write transmission conditions of this form on the outer 
boundary of Oh- The continuity condition (|45gJ) may be rewritten : 

W e o<S> l - W m \ v=1 = h 2 B m \ v=1 -v\ v=1 -v'o*!, 

where v and v e are respectively defined by (|39l) and (|37|) . We choose B m so that 
the right-hand side of the above equality vanishes: 



(48) 



v\ v =i -« e o$i _ 

^2 U 2 1*7=1) 



and thanks to Lemma [6] estimate ([42]) , the right-hand side of ([48]) is bounded in 
H 2 (T). 
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The condition (I45fl) is rewritten into: 



(49) 



— cW e o$i - -L d v W n \ =1 = — {d v uf + hd n B m \ v=l ) 



— 5„uo0$i - — 9„uf o$i, 

Me Me 



Observe that: 



d„< o $1 = 8 n u e o $ + hdlu e o $ + /i 2 / (1 ^ dlu e o $ dt, 

Jo * 

3 n uf o $! = 9„uf o% + h [ (1 - t)^ug o $ dt, 

Jo 

and recall that 

Me 

Therefore, in order for the right-hand side of (|49j) to be of order 2, we impose 
— d v B m \ v=1 - — ^«go$ - — 3„^o* = 0, 

Mm Me Me 

which implies 

(50) b m = rm - m + — d%v%o$ + — d n ulo$ , 

Me Me 



thanks to (j4"4"|) . We infer from (|4"5)) and (f5T))) that a m and b m are bounded in 
H 2 (T) independently of h, and therefore, since B m is polynomial in 77, it belongs to 
([0, 1]; # 2 (T)) , the space of functions, which are < ^°° in 77 € [0, 1] with values 
in H 2 (T). Particularly, there exists C > independent of h such that 

(51) V»/€[0,1], IIB"^, OH^-mt) < <7||0|| H . ( en). 

Observe that with such £? m , we have: 

a^B m | n=0 = 6 m + m , 

, Mm o2 e , Mm o e 1 

= mi H 9,>qO$o H <9„WiO$o, 

Me Me 



thanks to (|50p. Transmission condition (|28p with definition (|25|) of mi imply: 

Mc 



Therefore, transmission condition (I45c[) is simplified into: 

(52) —d n W C O$0 = J- 9^ m |^0 " ^9nS C O$ , 

Mc ^/^m 

It remains to define i? c . It is simply define by 

AB C + z c S c = 0, in O c , 
B e \ d0a = B m o$o 1 . 

Since B m 6 ff 2 (T), a classical argument and estimate ([5T|) imply that there exists 
C > independent on ft, such that: 

(53a) \\d n B c \ d0a \\HHdo a ) < C\\(f>\\ H >(dn), 

(53b) \\B c \\ H2+1/2{0a) <CU\\ H s (gn) . 
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To simplify our notations, we define 

g( V ,e) = ±;J2>(uZ + hu? + h 2 B m ), 

9^°) = i (— ( d n u ? + hd n B m \rj = i) - — d n u e Q o^! - — S„u?o$i 

ft \ Mm Me Me 

We equip L 2 (C) with the ordinary norm 

\H\l H c)= U 1 \a\ 2 d6dri) ' , 
and L 2 (T) with the ordinary norm 

II7Hl*(t) = ^h\ 2 de) 1 . 

We have chosen B m and B c such that there exists an ^.-independent constant C > 
such that 

(54a) IblU^c) < C\\(j)\\H*{dQ,), 

and 

(54b) ||5ilU 2 (T) < C\\(j)\\ H s (dn) . 

We rewrite Problem (|45l as follows: 
(55a) &W e + z e W e = 0, in O e . h , 

(55b) AW C + z c W c = 0, in O c , 

V( v ,6) e [o,i] xT, 

(55c) JzW n = ~/i 2 5 , 

with transmission conditions 

(55d) — d n W c ot> a =-^— d v W m \ v=0 -h 2 d n B c o^ , 

(55e) W c o$ = W m \ rj=0 , 

(55f) — d n W e o^x =-J- S# m | i + liV, 

(55g) VF e o$i = M /m |^ =1 

and the boundary condition 

(55h) 3„PH af2 =0- 

Now we are ready to perform L 2 estimates as it has been performed in [15] . In O c 
parameterized by Euclidean coordinates, the L 2 norm of a 0-form j3, denoted by 
l|/3|U 2 (O c ) 5 is equal to: 

\\P\\a°l2(o c ) = ll/3|U 2 (o c ), 
and the L 2 norm of its exterior derivative d/3, denoted by ||du|| a 1 l 2 (o c ) is equal to 

||d/3||AiL2(o c ) = ||V/3||i2( 0c) . 
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In O e ,h parameterized by Euclidean coordinates, the L 2 norm of a 0-form 7, denoted 
by ||7lU 2 (e> e , h ), is e q ual to: 

||7l|A°i 2 (o e , h ) = ||7ll-L 2 (e> e ,h)> 

and the L 2 norm of its exterior derivative dv, denoted by ||du||AiL3(e> e h \ is equal 
to 

l|d7l!AiL2( 0c-h ) = ||V7||z,2(o e> ). 

We multiply equalities (|55a|) . (|55b|l and (|55c() respectively by the conjugates of W e , 
W c and W m . Using transmission conditions (]55d|> — ( f55gl ) , we integrate by parts 
and we take the imaginary part of the result. To simplify the notations, we define 
II^IIao^^) and ||dW||^ 1Z;2(n) as follows: 

\\W\\l° L * m = \\W e \\lo LHOe:h) + \\W m \\l 0LUc) + \\W c \\lo LHOc)l 
lldl^lli^^o) = \\dW e \\l 1LHOeh) + \\dW m \\iiLUc) + W dwC W^L H o c ), 
Defining 

a = min(3?(z e ),3?(z m ),3(z c )), 



we obtain: 



d n B c W dvolr 
(1 + hn) ff iW e o$id6» 

1 



°Wtw(n) < — 

h 2 
H 



gW d? ? d(9 

Therefore, there exists an /i-independent constant C such that: 
\\W\\ 2 A o mn) < Cb?' 2 (Vh(\\d n B c \\ L2{ro) + \\gi\\ L2m ) 



(56) 



+ (11 Iffl^d^ \{\\W\\ A o L 2 m + \\dW\\ A1L 2 m ) 



Observe that C depends on the dielectric parameters and on the geometry of the 
domains. One more time, we multiply equalities ()55a|) — i|55b[) — fl55c[) respectively by 
the conjugates of W e , W c and W m . Using transmission conditions we integrate by 
parts and we take the real part of the result. Defining 



a = max{3fj(z e ),3?(z m ),3?(z c )}, 



Wc infer: 



\\dW\\ 2 A1LHn) <a\\W\\ 



+ 



31 W o$id<2 



d n B c W dvolr 

1 i-2-ix 



gW dr)d6 
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Using ([56]) we infer that there exists an /i-independent constant C > such that: 
\\W\\ A o LHn) + ||dW|| A i La( n) < Ch 3 / 2 (Vh (\\d n B c \\ L2{ro) + \\gi\\ LH r)] 

\g\ 2 d v d0 



o Jo 

Therefore, estimating the right-hand side of (l57|) with estimates (|53| and (|54|) . we 
infer that there exists an h- independent constant C > such that: 

||W|| A OHX (n) <CTl 3 / 2 ||0|| HS(9o) . 

Since we have respectively in the cylinder C 

u m - ug* - /wif = W ro - /i 2 B m , 

inO c 

?i c — Uq — ^ = W c - h 2 B c , 

and in O e 

w e - - ftu| = 

we have proved Theorem [TJ □ 

Remark 8 (Neumann boundary condition imposed on the cell). Consider the 
domain f2/j defined by: 

tt h = O c U o h . 

Let 7 be in H s (dflh), s > 7/2, and we denote by g and g the following function 
defined on the torus: 

(58) VfleT, s(0)=7o$i(0), 

(59) VxedO, gW^go^fi). 
Lei u be the solution of the following problem: 

V. ( — Vti J + qu = 0, inilh, 



Then, we have the following theorem: 

Theorem 9. We remember that ho is defined in Theorem^ Let 7 be in H s (dtth) 
We denote by Uq, u™, u\, and u™ the functions defined as follows: 

{Auq + z c Uq = 0, in O c , 
d n Uo\r = (Mc/A*m)fl, on T . 

In the membrane, the field is equal to: 

y(rj, 6) G [0, 1] x T, U ™=<o$oW- 
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The field u\ is the solution of the following problem in O c : 
(Aul + z c u\ = 0, in O c , 

[(/•WMc) d n u c 1 \ ra = (£g + <9 t 2 w§| ro ) + z m u c \ To , on T . 

In the membrane, we have: 

V(*7, 0) £ [0, 1] x T, u[ n = ?77 + u\ o $ . 
Let W be the function defined on flh by: 

u — (uq + hu'i) , in O c , 
u - (u™ o + fiu? o ,inO h . 
Then, there exists an h-independent constant C > such that 

\\W\\m iQh) < Ch 3 / 2 \\ Q \\ HH90c) . 

Remark 10 (The case z m = in Oh)- In Theorem^ we can replace z m by zero. 
The proof is then very similar, except that we need the following ineguality. 

Proposition 11. Let h be as in Theorem^ Let u be a function of class C 1 ([0, 1] x 
T). In the cylinder [0, 1] xT, we use Euclidean metric ([7]) written in local coordinates 
defined at Section QJ that is 

h 2 dn + (1 + hnn)d0. 



Then, there exists an h-independent constant C such that 



(60) 



IMIaol^c) < c (u^IIa^c) + I °^ d0 ) ■ 



Proof. Actually, according to the definition of ho in (O there exists two constants 
C\ and C2 depending on the domain O such that the following inequalities hold: 

r-l r 2ir 



(61a) 
(61b) 



Ia°l^(c) 



< C x h 



Jo 



\u( v ,9)\ z dOdn, 



Let us denote by (u) k for k £ Z the fc th -Fourier coefficient (with respect to 0) of \ 

u k = f u{6) e - 2 ™ k l L 66. 
Jo 

= 2mkuk, it is easy to see that: 



Since 



Vfc ^ 0, 
Using the following equality 

we infer 



Uk(rj)\ drj < 47r 2 



dry. 



d n u) (s)ds + Qq(0), 



\u (r])\ dr?<2 / (^A^fa) + 2\u (0)\dn. 
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Wc deduce directly inequality ([60]) . □ 

6. Application to the biological cell 

In biological cells, the membrane is insulating (see Fear and Stuchly [IT] or 
Sebastian al. [17j). This means that at mid frequencies, the ratio \z m \/\z c \ and 
|2m|/|*e| are small compared to h. Actually, the thickness is of order 10" 3 , while 
l^ml/l^d is about 10~ 5 (see [E]). We say that we work at mid frequency since 
we suppose that z c and z e are of order 1. Moreover, the relative permeabiltiy is 
constant equal to 1, thus in the following, wc suppose: 

The following results show that the asymptotics obtained by replacing z m by 
zero in the expansions of Theorem [T] give a good approximation of the electric field 
in the biological cell. We have the following proposition. 

Proposition 12. Let z c and z e be complex constants with strictly negative imagi- 
nary part. 

We suppose that \z m \ tends to zero and that there exists a constant c > such 
that: 

(62) 0< -A^L <c . 



Let z and z be such that: 



z e , in &:h , ( z e , in O e>h , 

Zm, in Oh, z= < 0, in Oh, 
z c , in O c , [z c , in O c . 

Let 4> in H 1 ^ 2 (dfl). Let u the solution of the following problem: 
(63a) Am + zu — 0, in Q, 

(63b) d n u\gn = cj>, in fl, 

and let v be such that 

(63c) Av + zv — 0, in CI, 

(63d) d n v\oQ = tp, in fl. 

Then, there exists a constant C such that: 

\\ U - V \\h^{Q.) ^ C 'l 2: mlll ( / , ll//i/2(^). 

Proof. First, using hypothesis (|62|) . we prove by classical argument that there exists 
an h- independent constant C such that: 

IM|ffi(n) < C\\(f)\\ H i/i (dn) . 

Then, we just have to write the problem satisfies by u — v in local coordinates in 
the membrane. As usual, we multiply in by u — v and we integrate by parts. Then, 
according to inequality ([60)) the following inequality holds: 

\\u- v\\ H1{Q) < C\z m \\\u\\ H i {n) , 

which ends the proof of the theorem. □ 

Using Proposition [12] and Remark I10[ we infer the following theorem. 
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Theorem 13. Let ho be in (0, 1) such that 

1 



h < 



l|/«Hoo" 

Let h be in (0, ho). 

Let z c and z e be complex constants with strictly negative imaginary part. 
We suppose that \z m \ = o(h) and that there exists a constant c > such that: 

Let cj) be m H s (dfL), s > 7/2. 

We denote by u the solution of the following problem: 

{Au + zu = 0, in fl, 
d n u = (f), on 90. 

We define (uq,Uq), u™, {u\,u\), and u™ are defined as follows. 

• The Q th order terms. The electric fields Uq and Uq are solution of the fol- 
lowing problem in O e .h U O c : 



(64a) 



Aug + z e u% = 0, in O eih , 
Auq + z c Uq — 0, in O c , 



with transmission conditions 

(64b) u c \r =u e \r , 
(64c) d n u c Q \ To = d n Uo\ ra , 

with Neumann boundary condition: 

(64d) d n u e \ dn = cj>. 

In the membrane, the field it™ is egual to: 
(65) V(?7, 6) G [0, 1] x T, u™ = u c o<S>o(0)- 

• The first order terms. The fields u\ and u\ are solution of the following 
problem in O e ^h U O c : 

!Au\ + z e u\ = 0, in O e> h, 
Au\ + z c u\ = 0, in O c , 
d nu\\ m = 0, 

with the following transmission conditions 

(66b) d n ul\r - d„Wi|r = d?u c \r a + d^u e \ To + Ad n u^\r , 

(66c) uj|r -ui|r„ = 0. 

Recall that M. is the curvature of dO c defined by §5§ . In the membrane, we 
have: 

(67) Vfa, 9) G [0, 1] x T, uf = rjd n u c o# +«iO$ . 
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Let W be the function defined in f2 by: 

!u - (uq + huf) j in C e ,h, 
u - (u(j + hu\) , in O c , 
u - (u^o®- 1 + hufo^- 1 ) , in O h . 

Then, there exists an h-independent constant C > such that 
\\W\\ HH n) <c(/i 3 / 2 + |z m |) W\\ H . m . 
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